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C03V2O8 is an orthorhombic magnet in which S=3/2 magnetic moments reside on two crystallo- 
graphically inequivalent Co^"*" sites, which decorate a stacked, buckled version of the two dimensional 
kagome lattice, the stacked kagome staircase. The magnetic interactions between the Co^^ moments 
in this structure lead to a complex magnetic phase diagram at low temperature, wherein it exhibits 
a series of five transitions below 11 K that ultimately culminate in a simple ferromagnetic ground 
state below T~6.2 K. Here we report magnetization measurements on single and polycrystalline 
samples of (Coi_a:Mga,)3V208 for a;<0.23, as well as elastic and inelastic neutron scattering mea- 
surements on single crystals of magnetically dilute (Coi_a:Mg2,)3V208 for i=0.029 and a;=0.194, 
in which non-magnetic Mg'^"'" ions substitute for magnetic Co^"*". We find that a dilution of 2.9% 
leads to a suppression of the ferromagnetic transition temperature by ~15% while a dilution level of 
19.4% is sufficient to destroy ferromagnetic long-range order in this material down to a temperature 
of at least 1.5 K. The magnetic excitation spectrum is characterized by two spin-wave branches in 
the ordered phase for (Coi_a;Mga;)3V208 (a;=0.029), similar to that of the pure x=0 material, and 
by broad diffuse scattering at temperatures below 10 K in (Coi_i:Mga;)3V208 (a;=0.194). Such a 
strong dependence of the transition temperatures to long range order in the presence of quenched 
non-magnetic impurities is consistent with two-dimensional physics driving the transitions. We 
further provide a simple percolation model that semi-quantitatively explains the inability of this 
system to establish long-range magnetic order at the unusually-low dilution levels which we observe 
in our experiments. 



I. INTRODUCTION 

Geometrically frustrated materials that are based on 
magnetic moments which decorate a lattice of triangu- 
lar networks have been of great recent interest due to 
their intriguing low-temperature magnetic properties and 
a rich variety of exotic ground states such as spin glasses, 
spin ices, and spin liquids which they exhibilP21_ 

C03V2O8 belongs to the M3V2O8 family (M=Co, Ni, 
Cu, Mn^PlSl which incorporates into its structure an 
anisotropic variation of the ideal two dimensional (2D) 
kagome lattice. The kagome lattice is itself a network of 
corner-sharing triangles. Within this structure, the mag- 
netic moments are carried by 3d transition metal M^"*" 
ions that decorate kagome layers. These layers buckle in 
and out of the plane, forming the kagome staircase struc- 
ture. C03V2O8 (CVO) crystallizes in the orthorhombic 
space group CmcaP. Edge-sharing Co^+Oe octahedra are 
situated slightly below and above the a-c plane, stacked 
along the crystallographic b-axis and are separated by 
nonmagnetic V^+04 tctrahedra. A representation of the 
structure including only the Co^^ ions is shown in Fig. 
la. The magnetic interactions between the Co^+ S=3/2 
moments are mediated by Co-O-Co superexchange, the 
resulting magnetic exchange pathways are indicated in 



Fig. la. as grey bonds. Within the buckled kagome 
layers that form the a-c plane (Fig. lb.) the S=3/2 
magnetic moments reside on two crystallographically in- 
equivalent Co^"*" sites: crystallographic (8e) spine (s) sites 
(blue) run in chains along the a-axis, and these chains are 
linked by (4a) cross-tie (c) sites (red). 




FIG. 1. The crystal structure of Co3V208- a) View of the 
kagome staircase considering only the Co'^"'' ions (red and 
blue) which are stacked along the 6-axis. b) View of the 
kagome plane projected on the a-c plane with the crystallo- 
graphically inequivalent cross-tie (red) and spine sites (blue). 
The magnetic exchange interactions discussed in the text are 
indicated. 
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The low-temperature magnetic phase diagrams of all 
known members of the kagome staircase family are com- 
plex and show considerable diversity in_their rich be- 
havior, despite their isostructural nature'^^^. For exam- 
ple, Ni3V208 displays 4 different magnetically ordered 
states below 10 K and has attracted much interest due 
to the discovery of multiferroic behavior in one of its 
incommensurate phases at low temperature"^ ^ ^ . The 
low-temperature phase diagram of C03V2O8 (CVO) has 
been studied by neutron diffract ion in both zercpKll and 
in finite applied magnetic field^^^. It has also been 
studied by optica l sp ectroscop^^^*^, by heat capacity 
and magnetizatiorPE21 as well as by /iSEp^ and by NMR 
technique j^^ * ^° l This extensive set of measurements has 
revealed a series of 5 different magnetic phase transitions 
below 11.2 K in zero applied magnetic field. These ulti- 
mately culminate in a simple ferromagnetic ground state 
below T~6.2 K. All five of the magnetic states display 
a preferred direction of the spins parallel to the a-axis, 
the easy axis of this system. From previous inelastic neu- 
tron scattering measurements, the exchange parameters 
of the microscopic Hamiltonian in the ground state of 
CVO have been extracted to some approximation, using 
linear spin wave theorjEl. These calculations approxi- 
mate the CVO system as two-dimensional and reveal that 
the exchange between magnetic moments on the cross-tie 
and spine sites is ferromagnetic and Jcs^l-25meV, while 
the exchange between adjacent spine sites (Jss in Fig. 
lb) is negligible. Significant uniaxial anisotropy terms 
in the Hamiltonian have also been found (~l-2 meV), 
consistent with the easy a-direction. 

Previous work has examined stacked kagome stair- 
case magnets in the presence of disorder, specifically 
(Coi_2,Ni3;)3V2 08, where random mixtures of Ni^+ and 
Co^"*" reside on the sit^^HH, The phase behavior of 
such systems is expected to be complex, as both Ni^+ and 
Co^"*" are magnetic, and a minimum of three exchange 
interactions would be necessary to describe the system, 
even if only a single M^"*" site was crystallographically 
relevant. 

Here we report on (Coi_j.Mg^)3V208 which can be 
thought of in terms of C03V2O8 in the presence of 
quenched magnetic vacancies, as Mg^'*" does not carry 
a magnetic moment. As the ultimate ground state in 
C03V2O8 is a simple uniaxial ferromagnet, one would 
expect this to represent an excellent model for a quasi- 
two dimensional Ising system in the presence of quenched 
disorder. As such, one would expect to be able to under- 
stand its phase diagram and ground state properties in 
some detail. 

Dilution studies of other Ising magnets have been car- 
ried out over a long period of tim^^. In most cases, these 
can be understood in the context of percolation theory. 
However, there are interesting examples where this is not 
the case. One such example is the case of the stacked tri- 
angular lattice Ising-like antiferromagnet in the presence 
of quenched non-magnetic impurities, CsCoi_j;Mgj,Br3, 
wherein the combination of antiferromagnetic exchange 



and triangular coordination leads the quenched magnetic 
vacancies to couple to the relevant order parameter as a 
random fielcP^. 

In this paper, we present magnetization measure- 
ments on both polycrystalline and single crystalline sam- 
ples of (Coi_2:Mga;)3V208 as wcll as elastic and inelas- 
tic neutron scattering measurements on single crystal 
(Coi_^Mg^)3V208 with X = 0.029 and x = 0.194. We 
study the phase diagram as a function of Mg^"*" doping, 
X, and show how a simple two-dimensional percolation 
model can provide an explanation for the drastic suppres- 
sion of the magnetic ordering at low temperatures in this 
material. In addition, we study the spin correlations and 
dynamics at low temperature as a function of doping, and 
show the existence of a Griffiths-like phasc!^ in the pres- 
ence of quenched disorder, when Tc{x) <T< Tc{x — 0). 



II. EXPERIMENTAL DETAILS 

Two large high-quality single crystals of 
(Coi_j,Mg^)3V208 with concentrations of a:=0. 029(3) 
and 0.194(4) were grown at McMaster University using 
an optical fioating-zone image furnac^^U^ Details of these 
growths and a subsequent X-ray structural refinement 
for these crystals are reported This work allows 

us to both quantify x, and to determine that the site 
dilution of the magnetic Co^"*" ions by non-magnetic 
Mg^+ is almost random. The single crystals resulting 
from these growths had a mass of ^ 8g and had approx- 
imately cylindrical shapes, with dimensions of 50 mm 
in length and 7 mm in diameter. Magnetization mea- 
surements were performed using a conventional SQUID 
magnetometer at McMaster University on several poly- 
crystalline samples with a range of Mg-concentrations 
{x — ~ 0.23) as well as on three single crystals {x=0, 
a;=0.029 and x=0.194). These single crystal samples 
were cut from the main crystal growths to a rectangular 
shape (~2.5 mm x 2 mm x 2 mm), and were aligned such 
that magnetization and susceptibility measurements 
could be performed with the magnetic field aligned along 
particular directions. All three single crystal samples 
used in these magnetization measurements had a mass 
of '^lO mg. Preliminary neutron scattering data was 
acquired using the C5 triple-axis spectrometer at the 
Canadian Neutron Beam Centre (CNBC) at Chalk River 
Laboratories. Further neutron scattering measurements, 
using both triple-axis and time-of-flight techniques, were 
performed at the NIST Center for Neutron Research 
(NCNR). The experiments were performed on two ^25 
mm long cylindrical (Coi_3,Mgj.)3V208 samples with 
x=0.029 and a:=0.194, which were aligned such that the 
[H,0,L] plane in reciprocal space was coincident with the 
horizontal scattering plane. 

Detailed elastic neutron scattering measurements were 
carried out using the cold triple-axis spectrometer SPINS 
at the NCNR, NIST. These measurements employed a 
vertically focussing PG-002 monochromator and fiat PC- 
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002 analyzer crystal with fixed final energy of Ef=5 meV. 
The coUimation used was [open, 80', 80', 80'] with an 
in-pile, cooled Be filter placed in the neutron beam in- 
cident to the sample in order to eliminate higher-order 
wavelength contamination. Inelastic neutron scattering 
measurements were carried out on the time-of-flight Disk 
Chopper Spectrometer (DCS) at the NCNR, NIST em- 
ploying a fixed incident wavelength of = 2.5 A, which 
allowed for the measurement of magnetic excitations up 
to energy transfers of AE ~ 10 meV. In this configura- 
tion of DCS, a resolution of ~ 0.9 meV was obtained at 
the elastic position. The samples were placed in a con- 
ventional ILL Orange cryostat with a base temperature 
of 1.5 K for both the SPINS and DCS experiments. 



III. MAGNETIZATION MEASUREMENTS 



The temperature-dependent mass susceptibility x(T) 
was measured in a field-cooled (FC) mode in an applied 
field of /ioH=0.005 T and this is shown in the inset of 
Fig. 2 for selected values of x. For the measurements on 
single crystals (open symbols), the crystals were aligned 
with their a-axis parallel to the applied magnetic field. 
Measurements on pure C03V2O8 show two well-defined 
magnetic transitions at Tn ~ 11-2 K and at Tc ~6.2 
K as indicated by arrows in the inset. We associate the 
higher temperature transition with a transition from the 
paramagnetic to an incommensurate antifcrromagnctic 
phase, and the lower temperature transition at ^6.2 K 
with the transition to the ferromagnetic ground state, 
as has been reported previousljff^. These two transition 
temperatures were extracted from x(T) ■ T curves, shown 
for a subset of the data in Fig. 3, resulting in the phase 
diagram shown in Fig. 2. 



We performed magnetization measurements using a 
SQUID magnetometer on a series of (Coi_a;Mga;)3V208 
samples with different doping concentrations in the range 
(a; = — 0.23). Measurements were carried out on poly- 
crystalline samples as well as on single crystal samples 
(x = 0, a; = 0.029 and x = 0.194). 




Mg concentration x in (Co,.^Mg^)3V20g 



FIG. 2. (Color online) The phase diagram of 

(Coi_iMga;)3V208 as a function of x obtained from 
magnetization measurements described in the text. The 
error bars result from the determination of the transition 
temperatures using data in Figure 3. The phase boundaries 
(dashed lines) are guides to the eye and are extrapolated to 
T = K. The theoretically expected site percolation limit 
for the perfect 2D kagome lattice {xc ~ 0.35)^^ is shown 
as the red star. The critical doping concentration x^ above 
which any magnetic ordering gets suppressed at T=0 corre- 
sponds to the determined percolation limit discussed later 
in the text and is shown as blue diamond. The inset shows 
representative scans of the field-cooled mass susceptibility 
X(T) for an applied field of ^oH=0.005 T, and the phase 
transition temperatures Tc and Tn for data points in the 
phase diagram. 
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FIG. 3. (Color online) Representative scans of the suscepti- 
bility x(T) • T reveal anomalies at the phase transition tem- 
peratures Tc and Tn (indicated by arrows). These transi- 
tion temperatures are used to construct the phase diagram 
shown in Fig. 2. Closed symbols refer to data obtained from 
polycrystalline samples, while open symbols represent mea- 
surements on single crystals. Note that the maximum in the 
X(T) ■ T curve in the x=0.194 doped sample is indicated as 
phase transition with a Tn ~ 5 K in Fig. 2, however, elas- 
tic neutron scattering data shows no clear evidence for long 
range magnetic correlations at any temperature above 1.5 K 
in this sample. 

These measurements reveal a suppression of both tran- 
sition temperatures with increasing x, and the complete 
suppression of any magnetic ordering at T=0 can be ex- 
trapolated to a critical doping concentration of Xc ^0.26. 
This result is somewhat unexpected since this doping 
level is lower than the theoretical percolation threshold 
for the destruction of long-range order for the ideal 2D 
kagome lattice, for w hich X c ~0.35 (indicated as red star 
in the phase diagram )PSl2ll. It is also much lower than the 
percolation threshold for any three dimensionally ordered 
system where Xc '^0.68-0.8CP^^2ll, Note that the perco- 



4 



lation threshold Pc reported in the hterature is related 
to Xc discussed here hy Xc = I — Pc- While Pc describes 
the critical concentration of magnetic ions below which 
magnetic long range order ceases to exist, Xc is used for 
the complementary description of the critical concentra- 
tion of non-magnetic ions above which long range order is 
impossible. We choose to speak of the percolation thresh- 
old Xc in this paper instead of pc as it is directly related 
to the concentration of non-magnetic Mg^+ ions in our 
doped CVO. 

From the phase diagram in Fig. 2, it can be seen that 
the ferromagnetic transition temperature Tc is very sen- 
sitive to magnetic dilution and changes almost linearly 
with X. It is lowered by ^^20% from Tc~6.2 K in the pure 
material to Tc^5.2 K for a doping level of a:=0.029, and 
for a;=0.194 it is suppressed below the lowest accessible 
temperature for our SQUID measurements (T^2 K). In 
contrast, the higher temperature Neel transition changes 
more gradually with dilution, exhibiting downwards cur- 
vature as a function of x. The reason for this difference is 
not clear, although we do note that while the transition 
at Tjv appears to be continuous, that at is clearly 
discontinuous even in the presence of quenched disorder, 
as we will discuss below. 



IV. ELASTIC NEUTRON SCATTERING 
RESULTS 



b) x~0.029 J _ 1.8K iio' 
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FIG. 4. (Color online) Contour plots showing the temperature 
evolution of the elastic magnetic scattering S(Q, E=0) around 
the Q=(0,0,2) Bragg peak position for a) a; = 0.029 and c) x = 
0.194. A background data set within the paramagnetic state 
at T = 20 K has been subtracted from the individual (0,0, L) 
scans so as to eliminate the weak nuclear component to the 
(0,0,2) Bragg scattering. Representative scans making up the 
contour plots are shown in panels b) and d). The instrumental 
resolution of SPINS is indicated in panels b) and d) as the 
black horizontal bar. The error bars represent one standard 
deviation in this and the following figures. 



Elastic neutron scattering measurements were per- 
formed on two doped single crystal samples (a;=0.029 and 
x=0.194) using the cold triple-axis spectrometer SPINS 
at the NCNR, NIST. The temperature dependence of 
the elastic magnetic scattering was followed around the 
(0,0,2) Bragg peak, which exhibits a weak nuclear compo- 
nent, and a strong ferromagnetic component to the Bragg 
peak below T(7^6.2 K in the pure, x—Q, material. The 
strength of the magnetic Bragg peak follows as a con- 
sequence of the form of the magnetic neutron scattering 
cross section which is sensitive to the component of mo- 
ment perpendicular to Q. Its strength indicates that the 
ordered moment points along a. A color contour plot con- 
structed from elastic longitudinal (0,0, L) scans through 
(0,0,2) as a function of temperature for both samples is 
presented in Fig. 4 a) and c). Representative scans are 
shown in Fig. 4 b) and d). Note that the intensity ap- 
propriate to the color contour plots in Fig. 4 a) and c) 
is on a log scale. 

In the data sets for both a::=0.029 and a::=0.194, a 
high temperature data set at T = 20 K (in the param- 
agnetic phase) has been subtracted, so as to eliminate 
the temperature-independent, nuclear scattering from 
the lower temperature data sets, and hence to isolate 
the magnetic scattering. It can easily be seen that the 
scattering profiles for the two samples are in striking con- 
trast to each other. As the temperature is lowered below 
T~Tc(a: = 0) ~ 6.2 K in the a;=0.029 Mg-doped sample 
(Fig. 4 a) and b)), diffuse scattering intensity charac- 



teristic of the development of short range correlations 
builds up, increasing in strength down to a temperature 
of Tc{x = 0.029) -5.2 K. Below -5.2 K it gives way to 
a dramatic increase in elastic scattering intensity which 
is now sharper in Q and characteristic of long range or- 
der. It is therefore associated with a phase transition to 
the ferromagnetic ground state. This can be seen more 
clearly in Fig. 5, which displays the temperature depen- 
dence of the net scattered intensity upon warming and 
cooling measured at the Bragg peak position (0,0,2) as- 
sociated with the order parameter, and at a Q position 
shghtly away from the Bragg peak, at (0,0,1.8), allowing 
for the parametrization of the diffuse scattering charac- 
teristic of the short range order. Based on the sharpness 
of the transition and the observed hysteresis in the order 
parameter, this transition is clearly of first-order nature. 

The diffuse critical scattering in the wings of the elas- 
tic Bragg peak is observed t o be drastically enhanced in 
a Griffiths-phase-like fashiorP^'^ in the temperature re- 
gion between Tc{x = 0.029) -5.2 K and T <Tcix = 
0) —6.2 K. Such a Griffiths phase shows enhanced spin 
correlations due to relatively rare, large correlated spin 
droplets'^'* that arise as a consequence of the random 
quenched disordei^^SI. The effect of thermal fluctuations is 
such that the system would be long range ordered in the 
absence of quenched disorder, and large percolating net- 
works of spins display enhanced correlations somewhat 
akin to long range order. 
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FIG. 5. The order parameter and diffuse scattering as a func- 
tion of temperature measured near tlie (0,0,2) magnetic Bragg 
peak in tlie a;=0.029 sample. Note tfie different intensity 
scales and counting times. 



The (0,0, L) scans of the magnetic scattering from this 
a;=0.029 sample are best described by a two-component 
hneshape consisting of a Lorentzian hneshape character- 
izing the short-range correlations and a resolution-limited 
Gaussian hneshape which describes the onset of long- 
range magnetic order. 

The results of this hneshape analysis are given in Fig. 
6, wherein the integrated intensities of the Gaussian 
and Lorentzian components of the magnetic scattering 
at (0,0,2) for the a;=0.029 sample are shown in panel a). 
A representative elastic scan with corresponding fits at T 
= 5.8 K, in the Griffiths-like phase, between Tc{x = 0) 
and Tc(x = 0.029) is displayed in the inset of panel a). 
For the higher doping level of the a;=0.194 sample, the 
scattering profile is in sharp contrast to that at a;=0.029, 
as shown in the color contour maps in Fig. 4 and with 
the data we now discuss in Fig. 6. In agreement with 
the magnetization data, no magnetic long range order is 
evident down to the lowest temperature measured, ~1.5 
K. Instead, the scattering profile is entirely dominated by 
a broad diffuse scattering component around the (0,0,2) 
Bragg position, indicative of very short-ranged correla- 
tions. After subtraction of the high-temperature nuclear 
Bragg intensity, it was found that a Lorentzian hneshape 
alone was appropriate to describe the magnetic scattering 
profile at all temperatures. Figure 6 b) shows the inte- 
grated intensity of the magnetic diffuse scattering around 
(0,0,2) and the inset shows several representative scans 
and the corresponding Lorentzian fits to this scattering 
for x=0.194. One can clearly see that the diffuse scat- 
tering is very weak at 10 K and that it develops notice- 
ably below Tc{x = 0) ~6.2 K, again consistent with a 
Griffiths-like phase between Tc{x = 0) and zero temper- 
ature, increasing monotonically on cooling down to 1.5 
K. 

The full widths at half maximum (FWHM) of the 
Lorentzian components characterizing the short-range 
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FIG. 6. (Color online) Elastic scattering in the a) a;=0.029 
and b) a;=0.194 doped (Coi_a:Mga;)3V208 samples, a) The 
integrated intensities of the Gaussian (black) and Lorentzian 
(blue) components to the Hneshape used to describe the mag- 
netic scattering in the a;=0.029 sample discussed in the text 
are shown. A representative scan with corresponding fits in 
the Griffiths-like phase at T = 5.8 K is shown in the inset, b) 
The integrated intensity of the Lorentzian hneshape used to 
describe the magnetic scattering in the a;=0.194 doped sam- 
ple is shown, while the inset displays representative scans in L 
along the (0, 0, 2) direction for different temperatures with fits 
to the Lorentzian hneshape as described in the text. For both 
panels a) and b) a T = 20 K high temperature non-magnetic 
background has been subtracted to isolate the magnetic scat- 
tering. Note that the peaks for the a;=0.194 doped sample in 
the inset of panel b) are much wider than the instrumental 
resolution at SPINS, which is shown as black horizontal bar. 



correlations in both the (Coi_a;Mgj:)3V208 samples with 
x==0.029 and a;=0.194 are compared in Fig. 7. It is clear 
that the FWHM of the x=0.029 sample drops to at 
Tc{x = 0.029) ~ 5.2 K where the sample develops long 
range order (as evidenced by a resolution-limited Gaus- 
sian hneshape), whereas magnetic scattering from the 
a;=0.194 sample maintains a finite width to the lowest 
temperatures measured, T — 1.5 K. At this tempera- 
ture, the (0,0,2) peak displays a FWHM of -0.25, cor- 
responding to a real-space spin-correlated-droplet within 
the kagome plane of about 24 A in diameter, approxi- 
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FIG. 7. (Color online) FWHM of the Lorentzian components 
shown for both doped single crystal samples. The dashed lines 
are guides to the eye. While for a;=0.029, the FWHM of the 
scattering near (0,0,2) drops to at 5.2 K as the sample de- 
velops long-range order, for the higher doping of a;=0.194, the 
scattering maintains a finite Q-width to the lowest tempera- 
tures measured. The trend for the evolution of the FWHM for 
a;=0.194 however, extrapolates approximately to at T=0. 



mately the size of two to three unit cells. While no phase 
transition to long-range order is observed for a;=0.194, 
the trend for the evolution of the FWHM of this mag- 
netic scattering approximately extrapolates to at T = 
0, indicating that x=0.194 is close to a quantum critical 
point. We also note that the FWHM for both samples ap- 
proach each other at the highest temperature measured 
~10 K, as we anticipate should occur. At these high tem- 
peratures, thermal fluctuations are expected to dominate 
the effects of the quenched disorder. 
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FIG. 8. (Color online) Maps of the elastic scattering inten- 
sity S(Q, E = 0) for (Coi_a;Mga;)3V208 with x=0.029 and 
x=0.194 and for two different temperatures, T — 1.5 K (top 
row) and T — 5.5 K (bottom row). 

We further investigated the nature of the elastic scat- 
tering using TOF neutron data taken on DCS at NCNR, 
NIST. The obtained maps of elastic scattering S(Q, E=0) 



FIG. 9. (Color online) Cuts through the color contour maps 
maps shown in Fig. 8 along 0,0, L. The observed (0,0,1) peak 
at T = 1.5 K in the a;=0.194 sample is a clear signature of 
incommensurate fluctuations of the form (0,5,L) that are ob- 
served in the x=0.029 sample at T = 5.5 K and T = 7 K. 
This is consistent with the phase diagram presented in Fig. 
2. 



are shown in Fig. 8 a), c) for a;=0.029 and in Fig. 8 b) 
and d) for x=0.194, each for two different temperatures, 
T — 1.5 K and T =5.5 K. One notes that the net elastic 
scattering at T = 1.5 K, following subtraction of a T = 20 
K data set, is different for the a::=0.194 sample than for ei- 
ther the pure a;=0 (not shown here) , or the lightly doped 
x=0.029 sample. Figure 9 shows that Bragg-like scatter- 
ing is observed for the a;=0.194 sample around all of the 
(0,0, L) Bragg positions for integer L, likely originating 
from short range incommensurate elastic magnetic scat- 
tering of the form (0, i5, L). This scattering, which does 
not lie in our horizontal scattering plane, is picked up 
in our DCS experiment due to the finite acceptance of 
the scattered neutrons out of the horizontal plane, and is 
consistent with the placement of the a:=0.194 single crys- 
tal sample on the phase diagram for (Coi_2:Mgj;)3V208 
as shown in Fig. 2. 



V. INELASTIC NEUTRON SCATTERING 
RESULTS 

We investigated the evolution of the dynamic spin cor- 
relations in (Coi_a;Mga;)3V208, the excitations out of the 
ground state, using inelastic neutron scattering measure- 
ments carried out on DCS at NCNR, NIST. Figure 10 
shows S(Q, E) observed in the three x=0, 2:=0.029, and 
x=0.194 single crystal samples along two perpendicular 
directions in reciprocal space. 

Panels a) and b) in the top row of Fig. 10 show S(Q,E) 
for the pure, x=0, sample at T = 200 mK, panels c) and 
d) correspond to the a;=0.029 sample at T = 1.5 K, and e) 
and f) correspond to a;=0.194. The respective 20 K, high 
temperature data sets for each x, have been subtracted 
for all data sets shown. The spin wave spectrum in all 
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FIG. 10. (Color online) S(Q,E) in x=0 (top), a;=0.029 (mid- 
dle) and a;=0.194 (bottom row) samples along two perpendic- 
ular directions in reciprocal space. All data shown has had 
its respective high temperature T = 20 K background data 
set subtracted. Note that the data has been taken at the base 
temperatures appropriate to each experiment, which was 200 
mK for x=0 and T = 1.5 K for both a;=0.029 and a:=0.194. 



three samples consists of two bands of spin excitations, 
a relatively narrow band near ~ 2 meV and a broader 
band between ^ 3 and 6.5 meV. These are consistent 
with earlier spin wave results on the pure x—0 sampl^Sl, 



• X 


= 


— ■— X 


= 0.029 


—1=1- X 


= 0.194 




-2 2 4 6 8 

Energy transfer (meV) 

FIG. 11. (Color online) Cuts through scattering inten- 
sity around [1.5,0,0] and integrating over a narrow range in 
II=[1.4,1.6] and L=[-0.2,0.2] on a logarithmic intensity scale. 
To facilitate the comparison between these two sets of mea- 
surements, the scattered intensity has been normalized to the 
incoherent elastic scattering of the pure sample. 

As expected from the elastic neutron scattering results 



discussed above, the inelastic spectrum for the pure and 
the x=0.029 samples look similar within their ordered 
states, although broadening of the spin wave excitations 
in energy is clear on doping, corresponding to shorter 
spin wave lifetimes in the presence of quenched impurites 
for the a:=0.029 sample. The inelastic scattering from the 
most heavily doped 2:=0.194 sample shows poorly defined 
bands of spin excitations which are most easily observable 
in the vicinity of Q=0. These are typical of the magnetic 
inelastic spectrum within the incommensurate phases or 
paramagnetic state at elevated temperatures. 

The inelastic magnetic spectrum can be further com- 
pared between the three x=0, a;=0.029, and x=0.194 
samples, by taking cuts through the color contour maps 
of S(Q,E) shown in Fig. 10, and normalizing to the in- 
coherent elastic scattering. This normalization should 
account for the somewhat different volumes of sample 
within the neutron beam. This is shown in Fig. 11, 
where we have shown an approximate constant-Q=(1.5, 
0, 0) scan, which has been constructed by integrating in H 
between 1.4 and 1.6 and in L between ± 0.2. The spectral 
weight of the sharp spin excitations falls off dramatically 
with increased concentration of quenched non-magnetic 
disorder. 
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FIG. 12. (Color online) The temperature evolution of the 
inelastic scattering S(Q,E) along [0,0,L], contrasting the 
x=0.029 and a;=0.194 samples. Panels a), b) and c) corre- 
spond to the a;=0.029 sample at the temperatures indicated 
in phase diagram in the top right panel. Panels e) and f ) cor- 
respond to the 2;=0.194 sample with temperatures as shown 
in the phase diagram. 

Further, we investigated the temperature evolution 
of the inelastic scattering in the two doped samples, 
x=0.029 and 0.194, as shown in Fig. 12 which shows 
S(Q,E) along the [0,0, L] direction in reciprocal space. 
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This is consistent with our elastic scattering studies, 
which placed the a;=0.194 doped sample at 1.5 K in a 
region of the phase diagram that is characterized by in- 
commensurate elastic magnetic order and rather diffuse 
spin wave excitations. This corresponds to the same re- 
gion of the phase diagram in which one finds the a;=0.029 
sample at elevated (T = 5.5 K and T = 7 K) tempera- 
tures. Note also the ressemblance of b) and c) with e), 
which are all placed in the incommensurate antiferromag- 
netic phase in the phase diagram of diluted CVO. 



VI. PERCOLATION CALCULATIONS AND 
DISCUSSION 

In light of the phase diagram obtained for 
(Coi_2,Mga;)3V208, shown in Fig. 2, and both the 
elastic and inelastic neutron scattering in the presence of 
quenched disorder, one can raise the question as to why 
dilution with quenched magnetic vacancies is so effective 
in destroying the magnetic order in this system. Figure 
2 shows the ferromagnetic transition to be suppressed 
to zero temperature at x ^0.19, while all vestiges of 
incommensurate magnetic order have been suppressed 
to zero temperature by x ^0.26. This is far below 
the percolation threshold for any three dimensional 
cooperative system, where critical concentrations of 
Xc ~ 0.7 {pc ~ 0.3) are typical^^— -. In and of itself, 
this shows that the weak three dimensional interactions 
along the third, stacking dimension, b, do not determine 
this criticality. These interactions are certainly weak 
compared to the strong interactions within the kagome- 
staircase plane, but they may also be frustrated as a 
consequence of the triangular coordination of spins on 
cross-tie sites which neighbor the spine-site spins along 
the b-axis. 

To further examine this question, we carried out a sim- 
ple two dimensional, zero temperature percolation calcu- 
lation relevant to (Coi_a;Mgj:)3V208. Our starting point 
for a calculation of the percolation limit in this system 
was a two-dimensional projection of the magnetic mo- 
ments within a single buckled kagome layer. For the 
calculations, a lattice of 100x100 unit cells was used, 
from which we randomly removed magnetic sites (ei- 
ther cross-tie or spine positions) with probability x. As 
Ramazanoglu et al.'J^ linear spin wave theory on the 
pure CVO material showed the magnetic interactions be- 
tween magnetic moments on the spine sites to be negli- 
gible, we then excluded magnetic interaction pathways 
between spine sites for the remaining steps in the cal- 
culation. We calculated the number of sites contained 
in the largest cluster, that is the number of sites con- 
tained in the largest near-neighbor connected region of 
the lattice, and checked whether any continuously con- 
nected path through the lattice (chosen arbitrarily from 
the left to the right edge of the lattice) could be found. 
If several such connected paths were found, the length of 
the shortest path through the lattice between sites mak- 



ing up the largest lattice-spanning cluster was calculated. 
We ran this simulation up to 30 times for each dilution 
probability x to get a good estimate of the percolation 
threshold, which we find to be at a;c^0.26. A calcula- 
tion for the full connectivity of the site-diluted lattice 
was performed as well, by including the interaction path- 
ways between spine sites. This calculation recovers the 
well-known percolation threshold for the two-dimensional 
kagome lattice of Xc^O.35 (pc=0.65). Though the fi- 
nite size of the system smears out the actual percolation 
threshold x^, the threshold for the original kagome lat- 
tice is well reproducecPSJ. Increasing the lattice size does 
not change the result for Xc significantly, and we there- 
fore conclude that our estimate of Xc ^ 0.26 for the case 
of (Coi_a;Mgj.)3V208, where the magnetic interactions 
between spine sites is neglible, is accurate. 



X (Mg concentration) 
B 0.0 0.1 0.2 0.3 0.4 0.5 0.6 




FIG. 13. (Color online) Percolation calculations showing the 
evolution of the largest cluster size in a lattice containing 
100x100 unit cells as a function of site dilution x. The in- 
set shows the length of the shortest path through the system 
connecting the edges of the lattice, if it exists. Above the 
percolation threshold, the lattice edges cannot be connected, 
and there is no shortest path. Note that the evolution of the 
largest cluster size is normalized to its value at zero site dilu- 
tion (a:=0). Error bars in the normalized largest cluster sizes 
are determined from the standard error. The determination 
of the shortest path through the system was based on only 
those data sets that showed connectivity of the lattice, and 
the error bars are given by the standard error of these data 
sets. 
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A summary of both simulations is shown in Fig. 13. 
In the top panel, the number of sites within the largest 
cluster, normalized to the total number of sites occupied 
by magnetic moments (which depend on x), is shown. 
The black squares denote the data for the case of full 
kagome connectivity, while data in blue circles represents 
the case in which spine-spine connectivity (ss-ss correla- 
tions) is excluded, corresponding to the physical picture 
relevant to (Coi_a;Mga;)3V208. In the inset, the length 
of the shortest path, if it exists, is shown for both cases. 
The point Xc for which the lattice cannot be connected 
and thus no shortest path can be found, is indicated by 
a dashed line and is determined as Xc ~ 0.35 for the full 
kagome connectivity and as Xc ^ 0.26 for the case ap- 
propriate to (Coi_j;Mg2;)3V208. It is interesting to note 
that the length of the shortest path increases relatively 
gradually for the full kagome connectivity, while exclud- 
ing spine-spine connections leads to a flat dependence 
on the length of the shortest path through the system, 
until a dilution level of about x ~ 0.20 is reached, be- 
yond which the length of the shorest path rises rapidly 
towards Xc ^ 0.26. This seems reasonable as the ran- 
dom removal of cross-tie and spine sites is equivalent in 
the fully connected kagome case, while the removal of 
cross-tie sites, which make up 1 /3 of all sites is more ef- 
fective than the removal of the majority spine sites, for 
the case where magnetic interactions between spine sites 
is negligible. This is because removal of a crosstie site 
"disables" a possible path to 4 spine sites, while removal 
on a spine site only "disables" the interaction path be- 
tween 2 crosstie sites. 

Representative lattice configurations at dilution levels 
of x=0.15, a;=0.30 and x=0.4:5 are shown in the bottom 
panel of Fig. 13, for both cases. The largest cluster in the 
lattice is colored in blue and the shortest lattice spanning 
path(s) are shown in red. As one can clearly see, at a dilu- 
tion level of X = 0.30, there exists a largest, lattice span- 
ning cluster in the case of full kagome connectivity (bot- 
tom middle panel) , while for the case of excluding spine- 
spine connectivity as is relevant for (Coi_2;Mg2,)3V208, 
only a finite-sized cluster exists, which does not span the 
full lattice and thus does not allow for long range mag- 
netic order. 



VII. CONCLUSIONS 

We have presented the low temperature phase diagram 
of the kagome staircase system CVO in the presence of 
quenched disorder, showing the transition temperatures 
T^v and Tc as function of non-magnetic Mg^"*" doping 
in (Coi_a;Mga;)3V208. Wc have found that the magnetic 
properties of this material at low temperatures are very 
susceptible to quenched disorder and that a doping level 
of a; ~ 0.26 is large enough to suppress any long range 



order to below 1.5 K. 

Based on the geometric arrangement of the magnetic 
Co^"*" moments within the kagome plane, a simple two- 
dimensional percolation model has been used to under- 
stand the effect of magnetic dilution on the possible ex- 
change paths. For that purpose, we considered a flat, 
two-dimensional kagome lattice with site dilution x, and 
calculated the percolation threshold. Based on previous 
inelastic neutron scattering on CVO, which found the 
spine-spine magnetic interactions to be negligible com- 
pared with the crosstie-spine interactions^^, we calcu- 
lated the percolation threshold of Xc ~ 0.26 for the case of 
crosstie-spine connectivity (excluding spine-spine bonds) . 
As a reference, we confirmed the percolation threshold for 
the kagome lattice as Xc ^ 0.35 {pc ^0.65). The calcu- 
lation of the percolation threshold excluding spine-spine 
connectivity is found to be in good agreement with the 
observed phase diagram that puts the percolation thresh- 
old at ~26% percent doping level. 

For two samples with doping levels of a;=0.029 and 
x=0.194, respectively, we have performed elastic and 
inelastic neutron scattering measurements. The elastic 
neutron scattering measurements could be understood in 
terms of Griffiths-like phase fluctuations appearing be- 
tween Tc{x = 0) and Tc{x), wherein enhanced short 
range order is observed. For the case of the lightly doped 
x=0.029 sample, such short range correlations collapse 
into a long range ordered ferromagnetic state via a first 
order phase transition at Tc^5.2 K. The more heavily 
doped x=0.194 sample displays only ferromagnetic short 
range order to the lowest temperatures measured, and 
shows the co-existence of incommensurate fluctuations 
at base temperature, consistent with the phase diagram 
for (Coi_a;Mgj;)3V208 dcscribcd above. 

Our inelastic neutron scattering measurements on the 
lightly-doped a;=0.029 sample show spin waves compara- 
ble to those observed in pure, x = 0, CVO, albeit with fi- 
nite energy widths, and hence spin wave lifetimes within 
the ground state. The equivalent magnetic excitation 
spectrum observed in the a;=0.194 sample shows mag- 
netic spectral weight in the same energy regime, but no 
well defined spin wave excitations at any wavevevctor. 
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